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The collaboration has measured a deviation from the standard model (SM) prediction in 
the like sign dimuon asymmetry in semileptonic 6 decay with a significance of 3.2 a. We discuss 
how minimal flavour violating (MFV) models with multiple scalar representations can lead to this 
deviation through tree level exchanges of new MFV scalars. We review how the two scalar doublet 
model can accommodate this result and discuss some of its phenomenology. Limits on electric dipole 
moments suggest that in this model the coupling of the charged scalar to the right handed it-type 
quarks is suppressed while its coupling to the d-type right handed quarks must be enhanced. We 
construct an extension of the MFV two scalar doublet model where this occurs naturally. 



INTRODUCTION 

The D0 collaboration has reported a 3.2 a deviation 
from the standard model (SM) prediction of the like sign 
dimuon asymmetry in semileptonic b decay [T]. This ob- 
servation joins past anomalous measurements of B s — > 
J/ip(f> and B~ -> decays that can be interpreted as 
a pattern of deviations consistent with new physics con- 
tributing a new phase in B s ^ mixing (for a recent global 
fit and discussion see [I]). 1 

If the explanation of the like sign dimuon asymme- 
try measurement and these correlated deviations are not 
statistical fluctuations, then new physics interpretations 
of this pattern are of interest. General operator analyses 
have been carried out [TUHT2] and indicate that operators 
induced by scalar exchange with unenhanced Yukawa 
couplings and order one parameters in the potential (i.e. 
order one Wilson coefficients) the mass scale suppressing 
the operators of interest is a few hundred GeV. 

Such a low mass scale is challenging to reconcile with 
known constraints from flavour physics unless minimal 
flavour violation (MFV) [131115] is imposed. New physics 
(NP) models with MFV have the quark flavor group 
SU{3) Ur x SU{3) Dr x SU(3) Ql only broken by the 
Yukawa couplings. However this scenario does allow new 
phases and so provides a framework for explaining the 
anomalies mentioned above without giving rise to flavor 
changing neutral current effects that are in conflict with 
experiment. 



1 The observed 2.6 cr deviation from the standard model (SM) ex- 
pectation [2] in the averaged measurements of B~ -t tc per- 
formed at Belle and Babar 3 -6 correlates correctly with a new 
physics (NP) contribution of a phase to with a sign consis- 
tent with the NP phase implied by the D0 dimuon measurement. 
Such a NP phase also correlates with the expectation of a shift in 
sin 2 /3 extracted from B B — > J/i/jtfi compared to the SM expecta- 
tion [3 [8] and extractions from measurements in B s — > J/ip K s . 
Such a consistent deviation is also observed, its statistical signif- 
icance is 2.1 cr. Also see |9| for a discussion on the evidence for a 
NP phase in B^ and B s meson mixing. 



In this paper we discuss scalar models with MFV that 
can explain these anomalies. We first review how tree 
level exchanges of a neutral complex scalar in a simple 
two scalar doublet model can lead to enhanced CP viola- 
tion in the B q meson system and discuss the phenomenol- 
ogy of this model. We then show that limits on electric 
dipole moments suggest that the coupling of the charged 
scalar to the right handed w-type quarks is suppressed 
while its coupling to the d-type right handed quarks must 
be enhanced to be consistent with the data. We construct 
an extension of the MFV two scalar doublet model where 
this occurs naturally 2 . 



SET UP 



We will utilize the recent fit of [2] to determine the new 
contribution to B q — B q mixing (here q = s,d). This fit is 
consistent in its conclusions with an earlier analysis [10] . 
The D0 result (a b SL ) and the SM prediction [5] (A b SL ) 
are given by 
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where the number of X b b — > fi + /i + Y events is given 
by N^~ + for example. The quoted a b SL is a combination 
of the the asymmetry in each B„, denoted a 



bq 
"SL- 



Each 

of these contributions to a b SL can be expressed in terms 
of the mass and width differences (Mi2,Ti2) of the B q 
meson eigenstates and the CP phase difference between 



2 Of course models with scalar doublets that are not supersym- 
metric suffer from the well know naturalness problem of keeping 
the doublets light compared to the Planck scale. 
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these quantities cf> q as 



where a and j3 are color indices and 
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Naively one can effect the SM prediction through modify- 
ing M±2 or Ti2 and both approaches have been explored 
in the literature. Modifying the decay width significantly 
p~6l IT7] as an explanation is problematic 3 and we will fo- 
cus on MFV NP explanations that involve a NP contri- 
bution (that includes a new CP violating phase) to Mf 2 . 

The effect of NP on B 8 and Bd mass mixing can be 
parametrized by two real parameters, h q > and a q by 
writing 



M q 
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NP 
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where the new physics contribution to the mass mixing is 
related to the standard model value of the mass mixing 
by, 



(Mf 2 ) SM h q e> 



(5) 



The models we discuss have h s = hd and a s — ad which is 
generally expected in NP models that obey MFV. 4 This 
scenario is argued to be a better fit to the current data 
then the SM in [2J, which is disfavoured with a p- value 
of 3.1cr. In this case, the best fit values are h q = 0.255 
and 2a q = 180° + 63.4°. The best fit magnitude of the 
correction h q is small but its phase is large. 

For simplicity in this paper we treat perturbative QCD 
in the leading logarithmic approximation and evaluate 
the needed matrix elements of four quark operators using 
the vacuum insertion approximation at the bottom mass 
scale. At the i-quark mass scale, in the SM, the effective 
Hamiltonian for B q — B q mixing is, 

= (F t * V tb fC SM (m t )btr q 2b^ q l (6) 



Here S(m^/M^) ~ 2.35 is a function of m?/Mj^ that 
results from integrating out the top quark and TU-bosons. 
Using 
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we have after running down to the 6-quark mass scale 
that, 

« 2 ) SM = {V: q V tb fC^{m t ) V fkm Bq . (9) 

Here rj ~ 0.84 is a QCD correction factor, C q M (m b ) — 
r)Cl M (m t ). 

The models for new physics we discuss generate the 
effective Hamiltonian at the top scale 5 



(10) 



Running down from m b operator mixing induces the anal- 
ogous operator with color indices rearranged. However, 
its coefficient is very small and we neglect it resulting in 
the relation, C NP (m b ) ~ ri'C NF (mt), where r)' ~ 1.45 
|20) . Again using the vacuum insertion approximation at 
the 6-quark mass scale we arrive at, 

(Mf 2 ) NP ~ (V t * V tb fC^{m t ) f-^) rif B m Bq . (11) 



Comparing with Eq. |5]) 
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MINIMAL TWO SCALAR DOUBLET MODEL 



3 The decay width can be removed in the relation between mea- 
sured quantities under the assumption of small CP violation in 
NP induced tree level decays of B q 1101 1181 and the anomalous 
measurements can still be fit to finding ~ 3 a evidence for a de- 
viation from the SM |2lllOI . Also, an explanation of the like sign 
dimuon asymmetry through a NP contribution to |r^ 2 | would not 
necessarily explain the anomalies in B s — > J/ty<j> and B — > r V 
with the correct correlation. 

4 It has been proven in |27| that new CP violating effects can be 
larger in B s than in B^ in nonlinear MFV. This observation 
has recently been explored in a general operator analysis 1121 
which showed that enhancements of CP violation in B s mixing 
over mixing by m s /m^ requires contributions in the MFV 
expansion out to forth order in both the up and down Yukawas 
for operators induced by scalar exchange. The results on the 
neutron EDM using naive dimensional analysis (NDA) on page 
4 disfavour order one down and up Yukawas with CP violating 
phases for these operators, so such contributions are expected to 
be very small. For alternative estimates of the relevant matrix 
element not using NDA see |19| . 



We now discuss how the minimal two scalar doublet 
model with MFV can have enhanced CP violation in B q 
mixing due to tree level exchange of neutral scalars. 6 We 
denote by H the doublet that gets a vacuum expecta- 
tion value and by S the doublet that does not. The 
Lagrangian in the Yukawa sector is 

£y = u i R 9 u i i QL S H + $ R g D s i Q Lj rf (13) 
+tf R Y u j i Q Lj S + di R Y D {Q Lj S^+h.c. 



5 For QCD running we don't distinguish between the top scale, 
weak scale and the mass scale of the new scalars we shall add. 

6 Previous analyses focused on scalar exchange to explain the like 
sign dimuon asymmetry include |21H23) . Also see |24l [25] for 
some phenomcnological studies of models of this form. 
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where flavour indicies i,j are shown and color and 
SU(2)l indices have been suppressed. MFV asserts that 
any NP also has the quark flavour symmetry group only 
broken by insertions proportional to Yukawa matrices so 
that Y v ^Yj-)^ are proportional to 9jji- *-*ne can 
construct allowed NP terms by treating the Yukawa ma- 
trices as spurion fields that transform under flavour ro- 
tations as, 



9u V v g v V c 



9d -> V D g D vL 



(14) 



where Vjj is an element of SU(3)u R , Vd is an element 
of SU(3)d R : and Vq is an element of SU(3)q l , i.e., the 
Yukawa matrices transform as gu ~ (3, 1,3) and go ~ 
(1, 3, 3) under the flavour group. MFV can be formulated 
up to linear order in top Yukawa insertions, or extended 
to a nonlinear representation of the symmetry |26| 12 7j . 
For enhanced CP violation in B q mixing we are interested 
in a nonlinear realization of MFV. It is sufficient to only 
expand to next order in insertions of gu so that 



Y 3 



Y D 3 i = VD9 D 3 i 



Vu9u t +Vu9u 3 k [(9u) i i9u)i] 
■v' D 9 D j M) k i(9u) l i 



We decompose the second scalar doublet as 



S+ 

s° 



(15) 



(16) 



where S° = (S R + iS$)/>/2. 
The scalar potential is 



V 



v 
~2 



m\ {S^Si), 



(17) 



+ (mlH^Si + h.c.) + Ai {H^Hi) (S^Sj), 
+ A 2 (fft< H 5 ) {S^ Si) + [X 3 H^ fftJ S . Si + h.c] , 
+ [X 4 H^ S^ Si Sj + A 5 S fl fftj ^ Hj + h.c] , 
+ A 6 (5 t4 5,) 2 . 

where i,j are SU(2) indices. Here v ~ 246GeV is the 
vacuum expectation value (vev) of the Higgs. Since we 
adopted the convention that the doublet S does not get 
a vev the parameters m\ and A5 are related by, 



0. 



(18) 



The spectrum of neutral real scalars consists of the Higgs 
scalar h and S R and Sj. However, these are not mass 
eigenstates. In the (h, S R , Sj) basis, the neutral mass 
squared matrix J\4 2 , where A3 is chosen real and positive, 



Where m% = m\ + (Ai + X2)v 2 /2. Within the convention 
that A3 is real, the couplings r)u, T)'u, t)d and rf D and A5 = 
X R + iX§ are in general complex. The mass eigenstate 
scalars N® with mass rtij are related to h, S R , Sj by the 
orthogonal transformations 

h = J20 hj N j: 
S R 



j 

E°«^> 



E°«^- 



(20) 



We find the CP violating NP contribution to B q — B q 
mixing from neutral scalar exchange is 



C NP (m t ) = (V2r]' D m b /v 
where F(x) = x 2 + ... , and 
A 



2 A 
2" 



(21) 



(22) 



For the rest of this paper we truncate the expansions in 
\f2m t /v at the leading non trivial term. So, for example, 
in Eq. (23 I we use F(x) = x 2 . 

For simplicity we now focus on the case where A5 = 0. 
Then h, S R and Sj are mass eigenstates and the term in 
the potential proportional to A3 is of interest as it leads 
to the mass splitting between the real neutral fields given 
by m R ~m 2 — 2 A3 v 2 . When this term is non vanishing 7 , 
the tree level exchange of Srjx generates a CP violating 
NP contribution to B q — B q mixing and 



C NP (m t ) = {n' D f 



V2 



m t 



\3 ml 



(23) 



In the above equation, the bottom quark mass mt ~ 
2.93 GeV is evaluated at the top quark mass scale and 



m S R/I =m s ± X 3 v . 



Using Eq. (12 1 the mass scale of the new scalars is 
given by 



^4 _ 20tt 2 A 3 |^|Vmg m| 
ms ~ h qV M^S{m? t /ml v ) 



+ AV 



Using the best fit value h q — 0.255 [2] we find that 
m% ~ (154GeV) 4 |^| 2 A 3 + (246GeV) 4 A2. 



(24) 



(25) 



Then for example with a value \rf D \ = 5 and A3 = 1 
the scalar mass scale is ms — 360 GeV. As the mass 



m 2 X^V 2 X^V 2 \ 7 Note that imposing custodial symmetry on the potential does 

• 1 9 ^ 9 1 no * f° rce ^3 ~~ ^ 0- Custodial symmetry violation is a measure of 

M = [ A5 U Trig + A3U . (19) the total mass splitting (m| - ) (m? _ ^ ) x (A| - (2A 3 ) 2 ) i> 4 

AgW 2 TTlg — A| / in terms of the potential given in Eq. jT7l. 



4 



Mass of S R and S/ 



600 
500 
400 



^ 300 



200 
100 





'to 



FIG. 1: Mass of Sr (solid) and Si (dashed) as a function of 
ri' D for fixed A3. The upper (blue curves) are for A3 = 1 the 
lower (black) curves are for A3 = 0.5. 



splitting is significant we show in Fig. 1 the masses of 
the neutral scalars Sr, Si as a function of r/' D . Moderate 
enhancements of r{ D avoid a light neutral state. 

We have checked that the mass scale ms required for 
B q — B 9 mixing is compatible with the constraints from 
K — K mixing. This compatibility is due to MFV, which 
causes the ratio of the relevant Wilson coefficients to scale 
asm 2 /m 2 \V td /V tb \ 2 . 

Next we derive constraints on r\\jr\D that come from 
one loop Feynman diagrams with charged S scalar ex- 
change. We will show that limits on electric dipole mo- 
ments imply that | Im [rjjj rjo ] | < 10 . Note that writing 
this as a constraint just on r\ur]r) depends on truncating 
a function of \f2m t jv at leading order. We also examine 
the constraint on Re^fy^] coming from experimental 
data on weak radiative B decay. 



Neutron Electric Dipole Moment 

The large CP violating phases needed in this two scalar 
doublet model contribute to other CP violating observ- 
ables. Notable among them are electric dipole moments 
(EDM's). We will restrict our discussion here to the dom- 
inant contribution that is not suppressed by small quark 
masses when naive dimensional analysis (NDA) [28 is 
used. It comes about through the colour electric dipole 
moment of the b quark 130] due to the effective Hamil- 
tonian 

SH bg = C gb g 3 m b ba^T a G a Xc7 e^^b, (26) 
inducing the dimension six CP violating operator 

Og — 93 fap-y^ VX ° Ga^pGp^GyXa, (27) 

of Weinberg [31] when the b quark is integrated out. Our 
discussion will largely parallel the discussion of [52 ■ As S 



couples both to the up and down type quarks it induces a 
one loop contribution to the effective Hamiltonian above 
with 



C gb (m s ) 



-lm[i]1j r/* D ] ( V2m t 



64:ir 2 m 2 s± 



f{m 2 Jm 2 s± ), (28) 



where 



f(x) = 



x — 3 



log a: 

(x-l) 3 ^2(x-l) 2 ' 



(29) 



Running to fi ~ m b using |29[ 130] and estimating the 
matrix element of the operator with NDA 8 gives in e-cm 
units 



d n ~ 2Im[?7* V* D }f(m 2 /m 2 s± ) 



ITeV 

m s ± 



10 



-26 



(30) 



This is a significantly larger effect on EDM's than 
quoted in the general operator analysis [ITj examining 
the effects of four Fermi operators on EDM's as this con- 
tribution is not suppressed by small mixing angles or light 
quark masses. For m s ± = 360 GeV the neutron EDM 
experimental bound of d n < 2.9 x 10 -26 e-cm implies that 
|Im[77y^]| < 0.26 We plot the allowed |Im[77^7?|,]| as a 
function of mass for this NDA estimate in Fig. 2. 



Allowed By neutron EDM using NDA 
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FIG. 2: Allowed |Im[7jy»jJ,]| as a function of the charged scalar 
mass. 

This suggests that llmfry^Ty^]! (and the sum of the ef- 
fect of all other cross terms such as |Im[7/^-*77^*]| etc.) 
is also small. However, given the uncertainties from 
hadronic matrix elements and given the fact that the pa- 
rameters that enter the contribution to EDM's are not 
identical to those in B q — B q mixing it is difficult to draw 
precise conclusions on the parameters in the model that 
are important for mixing. 



We use method (a) of [32] with a s (fi = I GeV) ~ 4 jr. 
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B — > X s 7 Constraints 

Of course the two scalar doublet model also gives new 
contributions to quantities that are not CP violating. 
Here we briefly review the constraints on this model from 
B — > X s 7 with these assumptions. The extra term in the 
effective Hamiltonian arises from charged scalar exchange 
and has the form, 



<5 W-B^tXs 7 



[V; s V tb ]C, (^s L a^F^b R ), (31) 



where e < is the electric charge. The Wilson coefficient 
is given by 



— Vu Vd 



2m 2 \ / 7 (m?/m|±) 



3 m| ± 



with 



1 (I 

4 



2 x \ogx 



(1 -x) s 



1 + logs; 
(1-z) 



(32) 



•(33) 



This operator's contribution to the measured branching 
fraction BR(B — > X s j)e„ >i.6GeV is known [33] 



BR(B^X s -f) E ^ 



1.6 GoV 



10" 



= 3.15 ± 0.23 -4.0u 2 Cv 



This constraint includes the effect of running this op- 
erator down to the scale mi,. Comparing to the world 
experimental average |34) we obtain a ler bound on the 
parameters of the form 



0.17<i?e[77^]/ 7 K/m 2 s± ) 



3 m 



, < 0.07. (34) 



For m s ± = 360 GeV we find -1.7 < Refaj^] < 0.7. 
Although this constraint is weak it is interesting that 
EDM's constrain lm[rj^ rf D \ while B — » X s j constrains 
Re[i^ rf D ). 



Collider Physics: Two Scalar Doublet Model 

Pairs of S particles can be produced through the tree 
level exchange of vector bosons produced through q q ini- 
tial states in the case of the Tevatron and LHC and e + e~ 
in the case of LEPII. 

From LEPII a bound on the mass scale of the new 
scalar doublet is obtained as no anomalous two and four 
jet events were seen when operating at y/s = 209 GeV 
where 0.1 of integrated luminosity was collected. 
The relevant cross sections in this case are given in |35j 
and the masses are bound to be 



m s ± > 105 GeV, m s o 



=o > 209 GeV. 



(35) 



We plot the allowed mg, A3 that satisfy this second bound 
for the minimal two scalar doublet model in Fig. 3. We 



have also performed an electroweak precision data fit. 
For scalar masses ~ 100 GeV the constraints are weak. 
The allowed mass splitting in this model is \mi — m± \ < 
200 GeV using the 95%CL region. 



Allowed By LEPII 
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FIG. 3: The allowed range of parameters ms, A3 in the two 
scalar doublet model considering LEPII direct production 
bounds. 

A light mass of Sj is allowed as these states must be 
produced in pairs through vector boson exchange and 
Sft can be heavy. However, a single neutral scalar parti- 
cle can be produced at the Tevatron in association with 
a charged scalar though W ± exchange. The partonic 
production cross section for producing S ± Sj or S ± 
(when the width is neglected) is 



G s 



1 - 



A/, 



w 



(36) 



where p is the center of mass momentum of one of the 
produced particles and s is the partonic center of mass en- 
ergy squared. We scan over the parameter space allowed 
by LEPII using this formula for the Tevatron produc- 
tion cross section (with MSTW 2008 PDF's [35]) where 
the W ± is produced off the valence u, d quarks. The 
renormalization scale in what follows is always varied be- 
tween mg/2 and 2 ms. The cross sections as a function of 
mass are shown in Fig. 4. The Tevatron can potentially 
constraint some of the allowed parameter space. Search 
strategies for pair production through weak boson fusion 
of charged S particles that decay into ibbt are also 
somewhat promising. In this case the production cross 
section for iris ~ 200 GeV is a ~ 1 fb with a signal of 
two b jets and two t jets is produced in association with 
tagging light quark jets at large px- 

At the LHC, production through the tree level ex- 
change of a vector boson is no longer dominated by W ± 
exchange. The cross sections for the pair production of 
scalars are all similar in their dependence on A3 and as 
a function of ms- We show <r(pp —> W ± — > S ± S R /j) for 
yft= 7TeV in Fig. 5. 
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Cross Section at TeVatron 
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FIG. 4: The cross section a(S ± Si) + a(S ± Sr) as a func- 
tion of m s for A 3 = (0.1,0.2,0.35,0.5,0.75,1) going left to 
right. Here we have also imposed custodial symmetry on the 
potential A2 = ±2A 3 for simplicity in the parameter scans. 

Cross Section at LHC 7 TeV 
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FIG. 5: The cross section a{S Si) + a(S ± Sr) as a function 
of m s for A 3 = (0.1,0.2,0.35,0.5,0.75,1,1.25) going left to 
right. Here we have also imposed custodial symmetry on the 
potential A2 = ±2A3 for simplicity in the parameter scans. 
The other production cross sections through Z* , 7* are simi- 
lar. 

Although these vector boson exchange cross sections 
for LHC are small, potentially observable signals at LHC 
do exist when t]d is larger than one, and Sr/i is made 
with large logarithms associated with collinear gluon 
splitting [37] and small pr of the spectator b quarks. 
The cross-section for the production of the lightest state 
b b — > Sj at leading log takes the form [35] 

6S \ u / J 1 x xs 

where b(x, (a) and b{x, fi) are the b quark and anti- 
quark PDFs respectively. The large logs from collinear 
gluon splitting are summed into the parton distribu- 
tion functions by choosing (i ~ mj. When we let 
T]d = \/ms/(154 GeV) and choose A3 = 1 the production 
cross sections for the LHC are given by Fig. 6. This pro- 



duction mechanism must compete with the large b pro- 
duction background from QCD. However, we note that 
this signal has a distinct feature in its reconstruction of 
a resonance in the highest px b quark pair with a larger- 
percentage of its total number of events at high px and 
small rapidity than the SM background, which has an 
approximate Rutherford scattering angular dependence 
in its production of b quarks. 



Sjt or S/ Production via b b 
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FIG. 6: The cross section a(pp — > bbSi) for Si produced 
through collinear gluon splitting and b quark fusion. Shown 
are the cross sections for y/s = 7, 10, 14 TeV. 



A MODEL WITH r)u NATURALLY SMALL 

The charged scalar in the two Higgs doublet model 
has couplings to the quarks that (at leading order in the 
Yukawa matrices) are given by, 

^charged = VuURgudLS + + ri D d R g D u L S~ + h.c. (38) 

We need large CP violating phases to get the fit value of 
a q . For large phases the limits on electric dipole moments 
suggest that \rjDr]u\ 10 _1 for charged scalars with mass 
of a few hundred GeV. Unless the charged scalars are 
considerably heavier than the neutral ones this bound 
is expected to hold in the minimal two scalar doublet 
model if the model is to give the central value for h q . If 
the limits on the electric dipole moments improve then 
this may become a more serious constraint. 

There does not seem to be any acceptable symmetry 
reason that this product is small and at the same time 
rjD is enhanced. In order to see that this is the case it 
is convenient to work in the basis where both H and S 
get a vevs vr and vs respectively. These can be chosen 
to be real. Then the charged scalar P + is the linear 
combination of the fields h + and S + , 



\M + V S 
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Because H no longer plays a special role we write the 
couplings of the scalars to the quarks as, 

C Y = e H uguQ L H + e s u R guQ L S (40) 
+ e' H d R g D Q L H^ + e' s d R g D Q L S^ + h.c. 

Here we are using MFV and taking the quantities that 
break the flavor symmetry to be gxjm ■ These matrices 
are proportional to the usual Yukawa matrices, 

Writing the charged scalar interaction as 



-'charged 



r)uUR9udLP + + VDd R gDU L P + h.c. (42) 



we find that, 



Vu 



VD 



e H v s - e s v H 



(43) 



One way to get r\u small while ?7_d is large is to have 
v R vs so that rju ~ —cs/ch and r\u ~ ^'h/ € 's an d 
take the corresponding ratios of e's to be small and large 
respectively. For their product to be small we also need, 
ese' H /e' s eH to be small. This is clearly possible, however 
there doesn't appear to be any symmetry reason behind 
these choices. 

Note that there is an interchange symmetry where 
H S that forces both r\u = r\D = 0. But in the 
limit of that symmetry, P + = (h + — S l+ )/\/2, and the 
symmetry's action on P + is P + — > —P + . Hence there is 
a stable charged scalar. 

The Glashow- Weinberg model [39] where H couples to 
the u-type quarks and S couples to the d-type quarks 
has, e' H = es = and so rju — vs/v R and rjjj = —vh/vs- 
In this model t/dVu = — 1> so when rju is small t\b is 
large, but their product cannot be made small even with 
a tuning of parameters. 

In this section we construct a MFV model that has rju 
small for a symmetry reason. New scalars that transform 
under flavour [3D] can naturally have a small r\jj. Con- 
sider a scalar field S% that transforms the same way as 
the Higgs doublet under the gauge group but as (1, 8, 1) 
under the flavor group, 



(44) 



We choose to represent the scalar in terms of the Gell- 
Mann matrices S 8 = S§ T a where a = 1 , . . . , 8 is a flavour 
index. The Yukawa couplings are given by 

£y = u R Yjj t (g^ (T a ) n {g D )\ Q Lj (45) 
+d R (T a r t (Y D y m Q L3 Sl a + h.c. 

where we have made the flavour indices explicit. We use 
hat superscripts to distinguish this model's parameters 



from the two scalar doublet model. Recall that, 
% \ = Vu 9u j i + Vu 9uM) k l (9u) l i\ + • 



Y D \ = f} D g D \ 



+ v' D 9 D \[(9 u ) k i(9u) l i\+----m 



The potential is given by 



2m\ Ti[Sl i S l 



8 ;J 



Ai H^ 1 Hi Tr[S« S Si ] + A 2 H 3 Tr[^' S Sl ] 
A 3 JTW Tr[sl S 8j ] + A 4 H^Tr^ S 8j S 8i ] 



A 5 H^Tr[Sl j S 8i S 8 A +h.c. 



(47) 



+ Ag Trf^ 1 5*81 S 8 j] + A7 Tr[Sg l S 8 j <5" S 8 i] 
+ AgTr^^si] Tr[4 J 'S 8j ] + A^Tr[S^ 8 ,]Tr[^S 84 ] 
+ X w Ti[S Si S 8j ] TrlSgSl'] + AuTr^^T^^] 

In the potential the index is an SU(2) index and the trace 
is over the down flavour index. We again rotate the phase 
of S 8 (relative to H) so that the A3 term is real, then the 
couplings and A^s and the rfs are in general complex. In 
the above potential there are no linear terms in S 8 after 
H gets its vacuum expectation value and so it is natural 
for it not to have a vev. 

In the potential and the Y's one can also insert ar- 
bitrary numbers of gp g D matrices between contractions 
of a down index. We work in the down basis so that 
go = diag(v / 2m- c (/v, \/2m s /v, y/2mt)/v). The interac- 
tions in the potential do not change flavour and are sup- 
pressed by ml/ v 2 so we neglect them. 

Keeping just the leading term in, \/2m t /v, the Wilson 
coefficient of the effective Hamiltonian as defined in Eq. 



(10) is 



C NP K) 



{f l ' D ) 2 {^2m t /v) i A 3 m 2 /6 



nil 



A§V/4 



(48) 



where m 2 s = 'in 
bound 



?+(Ar + A 2 ) t, 2 /4. 



This leads to the mass 



m z s ~ (98GeV) 4 |f^|'%+ (174GeV) 4 Ag\ (49) 

in terms of the parameters defined in the potential. The 
mass spectrum of the new doublet is given by 



m 2 s ± = m|- A 2 — , 



m 3 o = rh s + A 3 — , 



v 
~4 
v 2 
2 



1 - 2 \ u 

m s o = m s - A 3 y, 



(50) 



We show the masses of the neutral scalars for this model 
in Fig. 7. 
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FIG. 7: Mass of Sr (solid) and Si (dashed) as a function of 
rj' D for fixed A3. The upper (blue curves) are for A3 = 1 the 
lower (black) curves are for A3 = 0.5. 



This model has eight new scalar doublets. Never- 
theless, precision electroweak constraints are satisfied 
(when the Higgs is fixed to be m^ = 96^4 GcV for- 



ms 



> 



100 GeV) when |mj - m±| < 50GeV [35]. Con- 
versely, custodial SU(2) violation in such a light scalar 
doublet leading to a positive contribution to AT can raise 
the allowed mass of the Higgs in EWPD [551 HT] . 

In this model the coupling constant analogous to rju 
is naturally of order (m,b/v) 2 ~ 10~ 3 and the B — > Xs"f 
and neutron EDM effects of the model are suppressed as 
phenomenologically required due to MFV. 

The collider phenomenology in this model is very sim- 
ilar to the discussion on the two scalar doublet model. 
The LEPII constraints allow a larger parameter space 
due to the smaller mass splitting. The main differences 
for the Tevatron is that the cross sections we have dis- 
cussed are increased by an order of magnitude due to the 
larger flavour representation. Slightly smaller production 
cross sections through b quark fusion with low px spec- 
tator b quarks are expected at LHC as the normalization 
of the Gell Mann matrix decreases the cross section by a 
factor of three. 

We consider speculation on the UV origin of such a 
S$ doublet, or other accompanying non flavour singlet 
doublets that transform under the SU(3)u R as an 8 to 
be premature and beyond the scope of this work. 



CONCLUSIONS 

In this paper we discussed the new (i.e., beyond the 
minimal standard model) physics in the region of param- 
eter space for which the two scalar doublet model with 
MFV gives the additional contributions to B q — B q mix- 
ing that are hinted at by the data on flavor physics in the 
-B-sector. It requires additional light scalars that may be 
discovered at the Tevatron or LHC. Experimental limits 



on electric dipole moments suggest a region of parameter 
space that can occur naturally in some models where the 
new doublet of scalars transforms non-trivially under the 
flavour group. We constructed such a model. 

Acknowledgements 

We thank Z. Ligeti and Jim Cline for comments on 
the manuscript. This work was partially supported by 
funds from the Natural Sciences and Engineering Re- 
search Council (NSERC) of Canada. Research at the 
Perimeter Institute is supported by the Government of 
Canada through Industry Canada and by the Province of 
Ontario through the Ministry of Research & Innovation. 
MBW is greatful to Perimter Institute for their hospi- 
tality while this work was done. The work of M.B.W. 
was supported in part by the U.S. Department of Energy 
under contract No. DE-FG02-92ER40701. 



APPENDIX: EWPD CALCULATIONS 

The one loop results for the 5*8 model are the same as 
for the model discussed in [33 [35]. We use the STUVWX 
parameterization [33] of EWPD as the mass scale of the 
new scalars is ~ 100 GeV. The relevant results in terms 
of Passarino-Veltman functions [42 with standard defi- 
nitions 9 are 



6U ww (p 2 ) 



SU zz (p 2 ) 



<5IT 7 > 2 ) 



2ir 2 
1 



B 22 {p ,mj,m + ) + B 22 (p ,m R ,m + ), 

2\ 1 A C™2 \ 1 A t^2\ 



_sl_ 

2ir 2 c 



(l-2s 2 ) 2 (B 22 (p 2 ,m 2 + ,m 2 + ) 



1 



A (m 2 + ) 



B 22 (p 2 ,m 2 R ,m 2 j) - ^A (m 2 R ) - ^A (mj) 



2e 2 



1 



B 22 (p ,m + ,m + ) - -A (m + ) 



e 9l (l-2s 2 ) 



6U 7Z (p 2 ) 

For p 2 = these expressions become 



2 

2 ™2 



B 22 {p ,m + ,m + ) - -A (m + ) 



9l 



Sn ww (0) = ^ I i/(m +) m H ) + -f(m+,mi) 



SU zz (0) 
where 



8tt 2 c 2 



/(m fl ,m/) 



pi \ 2,2 2m{mi L m\ 

f(m 1 ,m 2 ) =m 1 +m 2 ^ 2 lo S ~~2 

mi — mi m. 



9 With c, s the cosine and sine of the weak mixing angle. 
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The derivatives of the vacuum polarizations are 



£n; 7 (o) = 

«i; z (o) = 
m' ww ( P 2 ) = 



B Q (0,m 2 + ,m 2 + ), 

B (0,m 2 + ,m 2 + ), 



6ir 2 

e.9i(l-2.s 2 ) 



2ir 2 



12tt 2 c 

" _ 1 A db 22 {p 2 ,m\,m\) 
6 <9p 2 

2 ™2 



2tt 2 c 2 



A- 

12 



dp 2 

db 22 (p 2 ,m 2 R ,m 2 j) 



2V- 



(l-2s 2 ) 



12 



-A 



dp 2 

db 2 2{p 2 ,m 2 +1 m 2 + 
dp 2 



Using these results we can construct the STUVWX 
parameters with the standard definitions [43] 



tS 



4s 2 c 2 



m zz (M 2 ) 



an zz (o) 



M 2 Z 



(c 2 



aT = 

all 
4^ 



s c 

<JIW(0) 



<5n' ZT (o) - <m 77 (o), 



<5iW(o) 



6Ilzz(M%) - m zz (o) 



Mf 



-s 2 OT' (0)-2scOT^„(0), 



alT 



«5n' zz (M 2 ) 



WW 



Ml 
JIWtM 2 ,) 



OT zz (0) 



Ml 



Here A is the divergence that cancels in the pseudo- 
observables STUVWX but we note we calculate in dimen- 
sional regularization and MS in d = 4 — 2e dimensions. 
As the number of degrees of freedom in this Sg, model and 
in the model [35] are the same, we can directly use the 
detailed fit results on the allowed masses (determined 
from these formulas) presented in [35]. These results 
generally allow masses for fixed rrih = 96±f 4 GeV when 
m s ?i 100 GeV characterized by \mj — m± \ < 50 GeV for 
Sg. 
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